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Summary 

The Earth Radiation Budget Experiment (ERBE) 
uses scanning radiometers to measure shortwave and 
longwave components of the Earth’s radiation field 
at about 40 km resolution. It is essential that these 
measurements be accurately located at the top of the 
Earth’s atmosphere so they can be properly inter- 
preted by users of the data. Before the launch of 
the ERBE instrument sets, a substantial emphasis 
was placed on understanding all factors which in- 
fluence the determination of measurement locations 
and properly modelling those factors in the data 
processing system. After the launch of ERBE in- 
struments on the Earth Radiation Budget Satellite 
and NOAA 9 spacecraft in 1984, a coastline detec- 
tion method was developed to assess the accuracy of 
the algorithms and data used in the location calcula- 
tions. Using inflight scanner data and the coastline 
detection technique, the measurement location errors 
are found to be smaller than the resolution of the 
scanner instruments. This accuracy is well within 
the required location knowledge for useful science 
analyses. 

Introduction 

The Earth Radiation Budget Experiment (ERBE) 
is a three-satellite mission designed to measure com- 
ponents of the Earth’s radiation field. Three satel- 
lites, Earth Radiation Budget Satellite (ERBS), 
NOAA 9, and NOAA 10, are now in orbit. Each 
satellite carries two ERBE instruments, the scan- 
ner and the nonscanner. The scanner has three 
Earth viewing radiometers, which make a total ra- 
diation measurement and broadband radiance mea- 
surements in the shortwave and longwave parts of 
the energy spectrum. The three scanner channels are 
coaligned and nominally scan the Earth every 4 sec- 
onds from horizon to horizon in a plane normal to the 
orbit plane. The nonscanner has four Earth viewing 
channels and a solar monitor channel. Two of the 
Earth viewing channels measure the total radiation 
and two measure the shortwave component. A single 
nonscanner measurement covers a large area of the 
Earth, and the scanner has a field of view (FOV) of 
about 40 km in diameter at nadir. 

Two different spacecraft are used to carry the 
ERBE instruments, and each is operated by a 
separate agency. The National Aeronautics and 
Space Administration (NASA) Goddard Space Flight 
Center (GSFC) operates the ERBS, which was 
launched on the Shuttle Challenger in October 1984. 
The Earth Radiation Budget Satellite is in a near- 
circular 600-km orbit with an inclination of 57° . The 
National Oceanographic and Atmospheric Adminis- 


tration (NOAA) controls the NOAA 9 operational 
weather-monitoring satellite, which was launched 
into a polar orbit in December 1984. The NOAA 9 
satellite is in a sun-synchronous near-circular 856-km 
orbit. Langley Research Center (LaRC) is responsi- 
ble for the ERBE instrument development, science 
management, and data processing. 

The development of the experiment and the sen- 
sors, as well as the stringent accuracy requirements 
of the mission, are described in reference 1. One 
element of the overall measurement accuracy is the 
requirement to locate the center of the FOV of each 
measurement. Both the surface characteristics at 
the measurement location and the satellite-Earth- 
Sun geometry are important factors used in the in- 
terpretation of the radiometric channels. Thus, it 
is essential that the measurement be accurately lo- 
cated with respect to the surface of the Earth. This 
task is complicated by the fact that, since there is 
no imaging system onboard the ERBS spacecraft, 
landmark identification is difficult. 

The ERBE project personnel began to analyze 
the factors affecting the accuracy several years be- 
fore launch of the satellites. Beginning with the 
source of the data, LaRC documented agreements 
with GSFC and NOAA which define the data in- 
terfaces to the LaRC processing system. The doc- 
uments specify a consistent set of orbital position 
data for each satellite and determine the coordinate 
systems used for measurements made from each 
spacecraft. The agreements also specify the format, 
content, and interpretation of the spacecraft attitude 
data derived from the onboard attitude determina- 
tion and control systems. The orbital position data 
and spacecraft attitude data permit the calculation 
of the location and orientation of the satellite, fac- 
tors essential to the accurate determination of the 
measurement location. Since the coordinate systems 
used to derive the various data types are different, 
position and attitude data must be transformed to 
a common coordinate system for interpretation and 
subsequent calculations. This step is the first key 
factor in understanding the accuracy of the measure- 
ment locations. 

Another accuracy factor involves the alignment 
of the radiometric instruments themselves with the 
satellite structure. Langley Research Center engi- 
neers monitored the instrument development con- 
tract with the prime contractor and measured the 
alignment of the instrument coordinate system with 
each spacecraft. These alignments are included in 
the measurement location calculations. 

Many of the remaining steps in calculating 
the measurement location require geometry 
transformations or modelling of the Earth’s surface 


to describe intersection points. All these steps were 
carefully derived and checked prior to launch. 

The purpose of this paper is to describe each 
step in the calculation of ERBE measurement lo- 
cations and to evaluate the accuracy of the cal- 
culations. Since the scanner instrument requires 
much more stringent pointing accuracy than the non- 
scanner instrument, the discussion is limited to the 
scanner instrument. The same concepts are used in 
locating the nonscanner measurements. The moti- 
vation for the analysis is to verify that the point- 
ing algorithms are correct and to confirm the va- 
lidity of the radiometric data location for science 
investigations. 

Symbols 


A,B,C,D,E,F ,G 

rotation matrices 

h 

height of top of atmo- 
sphere above Earth ellip- 
soid, km 

P 

vector from center of 
Earth to spacecraft 

R 

Earth radius, km 

V 

velocity, m/sec 

V 

vector from spacecraft to 
intersection with Earth 

V 

unit pointing vector 

X,Y,Z 

Cartesian axes 

a 

instrument azimuth 
angle, deg 

0 

instrument elevation 
angle, deg 

A 

difference between 
geodetic and geocentric 
latitude, deg 

e 

error, deg 

V 

spacecraft heading angle, 
deg 

d 

pitch angle, deg 

A 

longitude, deg 

V 

inertial longitude of 
Greenwich meridian, deg 

$ 

latitude, deg 

<i> 

roll angle, deg 


Subscripts: 

a 

C 

c 

D 

E 

e 

F 

G 

H 

/ 

N 

P 

P 

a 

0 

Abbreviations: 

ADACS 

ERBE 

ERBS 

FOV 

GSFC 

LaRC 

NASA 

NESDIS 

NOAA 

TOA 


yaw angle, deg 


along track 

geocentric 

cross track 

detector 

ERBS 

equatorial 

Earth fixed 

geodetic 

local horizon 

inertial 

NOAA 

pedestal 

polar 

azimuth 

elevation 


Attitude Determination 
and Control System 

Earth Radiation Budget 
Experiment 

Earth Radiation Budget 
Satellite 

field of view 

Goddard Space Flight 
Center 

Langley Research Center 

National Aeronautics and 
Space Administration 

National Environmental 
Satellite and Data Infor- 
mation Service 

National Oceano- 
graphic and Atmospheric 
Admininstration 

top of atmosphere 
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Calculation of Scanner Measurement 
Location 

The determination of the scanner-measurement 
location can be separated into two major steps. In 
the first step, the unit pointing vector of the optical 
axis of the detectors is calculated in a local-horizon 
system located at the spacecraft. In the second step, 
the detector pointing vector is transformed from the 
local-horizon system to an Earth fixed coordinate 
system in which the measurement location is deter- 
mined on a surface 30 km above the Earth. This 
surface is called the top of atmosphere (TOA) and 
has been selected as a reference surface for outgoing 
radiation. 

Calculation of Detector Pointing Vectors in 
Local-Horizon System 

In this first major step in determining the scanner 
measurement location, a unit vector which is aligned 
with the scanner optical axis is computed in a local- 
horizon coordinate system. This section describes 
the different coordinate transformations that relate 
the components of the detector pointing vector in the 
detector coordinates to its components in spacecraft 
coordinates. The discussion includes consideration of 
detector misalignments, a description of the teleme- 
try data, and the computations required to determine 
the pointing vector in spacecraft coordinates. This 
description is followed by a discussion of the transfor- 
mation of the pointing vector, through the use of the 
spacecraft attitude data, to the local-horizon coordi- 
nate system of the particular spacecraft. Finally, the 
pointing vector is transformed to a local-horizon sys- 
tem that is common to both the ERBS and NOAA 
spacecraft. 

Detector-to-Spacecraft Coordinate 
T ransfor mations 


the pedestal. The azimuth assembly (the entire sys- 
tem below the pedestal) can rotate (change azimuth) 
about the Xp- axis, so that the instrument can be 
pointed during solar calibration and so that mea- 
surements can be made while scanning at azimuth 
positions from 0° to 180°. The azimuth assembly is 
shown at the zero-azimuth (a = 0°) position, and 
a positive rotation (azimuth increasing) produces a 
positive rotation vector along the Xp- axis. The el- 
evation beam on which the detectors are mounted 
(small substructure on azimuth assembly) can rotate 
(change elevation or scan) about the Y^-axis. Ele- 
vation, or scan angle, is measured from the Y a -Z a 
plane, and a positive rotation produces a positive ro- 
tation vector along the Yq- axis. The detectors shown 
in figure 1 are at an elevation scan angle of 90°. 

To determine the complete matrix which trans- 
forms the detector pointing vectors into spacecraft 
coordinates, the matrix which transforms a pointing 
vector in one of the spacecraft coordinate systems 
to its components in the detector coordinate system 
should be defined first. The complete series of rota- 
tions which effect this vector transformation is given 
below: 


\X D 


'x E 


Y D 

Zd\ 

= DCBA# 

Y E 
. . 

(1) 


'X D ' 


X N 


Y d 

[z d J 

= DCBAjv 

Yn 
Zn . 

(2) 


where the pedestal-to-spacecraft alignment for each 
spacecraft is 


Figure 1 illustrates the coordinate systems which 
relate the position and orientation of the scanner 
detectors to the coordinate system of the instrument 
pedestal. The pedestal physically attaches the in- 
strument system to the spacecraft, and the relation- 
ship between the coordinate system of the pedestal 
and those of the TIROS N (NOAA) and ERBS space- 
craft is illustrated in figure 2. Only the orientation 
of these coordinate systems relative to each other is 
important to the measurement location calculations, 
and a unit vector describing the direction of the scan- 
ner detectors is assumed to have its origin at the focal 
point of the detectors. 

The three scanner detectors shown at the bot- 
tom of figure 1 have 2 degrees of freedom relative to 
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and where the azimuth, elevation, and detector rota- 
tions for the ERBS and NOAA spacecraft are 


r 1 
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0 

0 cos a 

sin a 

Lo - 

sin a 

cos a 

cos/3 

0 

-sin/3 

0 

1 

0 
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0 

cos/3 

■ 0 

0 1- 
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1 0 


.-1 

1 

O 
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The matrix A is the transformation from the 
spacecraft coordinate system to the instrument 
pedestal coordinate system. The alignment of these 
systems is shown for the ERBS and NOAA space- 
craft in figure 2. The matrix B is the transformation 
from the pedestal coordinate system to a new coor- 
dinate system (X Q ,Y a ,Z a ), which results when the 
instrument azimuth axes are rotated away from the 
pedestal axes about the Xp-axis. (See fig. 1.) This 
azimuth rotation is positive when the rotation is pos- 
itive about the Xp-axis. The matrix C is the trans- 
formation from the azimuth coordinate system to the 
coordinate system (X^F^Z^), which results when 
the instrument elevation beam is rotated away from 
the azimuth coordinate system about the F^-axis. 
This elevation rotation is positive when the rotation 
vector is positive along the F^-axis. In this new co- 
ordinate system, the detector scan head is aligned 
approximately with the Xp-axis when the elevation 
or scan angle is zero, whereas the elevation angle out- 
put of the instrument assumes that the detector scan 
head lies approximately in the Yp — Zp plane when 
the elevation angle is 0°. A final transformation, ma- 
trix D, transforms from the elevation coordinate sys- 
tem to the detector coordinate system by a “90° ro- 
tation about Yp to give the vector components in the 
detector coordinate system (Xp>, Yp), Zp>). 

Finally, to transform the pointing vectors from 
the detector coordinate system to one of the space- 
craft coordinate systems, the transpose of the matri- 
ces in equations (1) and (2) is used. The complete 
transformations of the vectors for the two spacecraft 
can thus be written as follows: 
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( 4 ) 


If there are no misalignments of the detectors with 
respect to the axes of the spacecraft on which they 
are mounted, the components of the pointing vector 
which is aligned with the detector optical axis, or 
initial pointing vector, are given by 

X D = 1 
Y d = 0 
Z D = 0 

Detector misalignments. This section describes the 
procedures that were followed during construction of 
the ERBS instrument to minimize and record detec- 
tor misalignment, and the methods by which these 
measurements were incorporated into the algorithms 
used to compute the pointing vectors. These pro- 
cedures were carried through the several contractors 
involved in the design and construction of the instru- 
ments and the spacecraft, and were carefully moni- 
tored by LaRC engineering personnel at all stages. 
Angular alignment data and analysis for the instru- 
ment on NOAA 9 have not been reduced to values 
which can be used to compute the initial pointing er- 
ror. Therefore, no detector misalignments have been 
assumed in computing the in-flight detector pointing 
vectors for NOAA 9. 

Preliminary alignment measurements for the in- 
strument on the ERBS spacecraft have been used to 
compute the initial pointing vector of the detectors. 
The components of that initial pointing vector are 

X D = 0.999999182 
Y d = -0.001197480 
Z D = -0.001005180 

More recent values of scanner detector angle mis- 
alignments for the instrument on the ERBS space- 
craft give an initial detector pointing vector with the 
following components: 

X D = 0.999999654 
Y d = -0.000824330 
Z D = 0.001221080 
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If the revised alignment data were assumed to be 
correct, the pointing vectors computed using the pre- 
liminary alignment data result in a 1.2-km scanner 
measurement location error at nadir in the direction 
of the scan motion. This error, which is a bias in 
the instrument and spacecraft coordinate systems, 
changes direction in the local-horizon system when 
the ERBS spacecraft yaws 180° approximately every 
36 days. The resulting along-track error using the 
preliminary alignment data is about 200 m down- 
range when the spacecraft is flying with the X-axis 
forward and 200 m uprange when the X-axis is rear- 
ward. Although there are small differences between 
the preliminary and revised values, the preliminary 
data have been used in all in-flight computations of 
ERBS scanner pointing vectors. 

Overview of telemetry data . The telemetry data 
from the ERBS and NOAA spacecraft are acquired 
and processed initially for LaRC by the Information 
Processing Divisions at GSFC and NOAA/NESDIS, 
respectively. Twenty-four-hour data sets are format- 
ted into logical data records, each containing 16 sec- 
onds of data and an associated universal time. The 
various types of data are sampled at different rates, 
and the universal time of any data word in a record 
is obtained by adding a specified time offset to the 
record time. The time offset for any data word in a 
record is a constant for all records for each spacecraft. 

The scanner instrument has a 4-second scan 
period, and each 16-second record contains 4 com- 
plete cycles of instrument radiometric and house- 
keeping data. The data in a 4-second cycle consist 
of 74 sets of measurements, and each set consists of 
one measurement from each of the three radiometric 
channels and a corresponding value for the detector 
scan head elevation angle. The instrument house- 
keeping data include, besides temperature and volt- 
age measurements, detector azimuth and elevation 
angles and instrument command and status informa- 
tion. Spacecraft housekeeping data include space- 
craft attitude data and spacecraft operational sta- 
tus information. The following sections describe how 
the raw instrument data and spacecraft housekeeping 
data are processed to determine the detector pointing 
vectors for the scanner radiometric measurements. 

Azimuth-angle data . In the normal Earth viewing 
data-acquisition mode, the scanner azimuth is set for 
crosstrack scanning (in a plane normal to the orbit 
plane). The azimuth position is 180° for the ERBS 
spacecraft and 0° for the NOAA spacecraft. Peri- 
odically, the azimuth assembly is rotated to obtain 
solar calibration data, and it is occasionally rotated 
to the 90° azimuth position to obtain along-track 
measurements. 


Four raw values (one for each scan cycle) of 
instrument azimuth position are contained in each 
16-second telemetry data record. Instrument com- 
mand and status data are checked to deduce the 
operational mode of the instrument and to determine 
whether a mode-change command has been issued for 
the record or whether the azimuth assembly is in mo- 
tion. If the operational mode is explicitly determined 
from the data, if no new operational command has 
been issued for the record, and if the azimuth assem- 
bly is not in motion, the first raw value of azimuth 
position in the record is decoded and converted to an 
azimuth angle. If the converted azimuth value is an 
acceptable azimuth angle and passes all edit tests, it 
is used in the computation of all space scanner detec- 
tor pointing vectors for the record. In all other cases, 
no valid detector pointing vectors are computed for 
the record. 

Elevation- or scan-angle data • The three radio- 
metric measurements and the associated elevation 
angle in a data set are all sampled within a few 
microseconds of each other. Sixty-two of the 74 mea- 
surement sets in a 4-second scan cycle are sampled 
every 33^ milliseconds. For the normal Earth scan 
mode, these 62 measurements are made while the de- 
tector is scanning the Earth from horizon to horizon. 
Even though each measurement set is sampled at a 
constant time offset relative to the record time, and 
even though the Earth viewing measurement sets are 
equally spaced in time, the elevation angle associated 
with a specific measurement set in a scan cycle differs 
from scan to scan within a record and from record to 
record. Because of this variation from scan to scan, 
a new detector pointing vector is computed for every 
measurement set. 

A detector pointing vector is computed only if the 
corresponding elevation angle is valid. The valida- 
tion of the elevation-angle data involves checking the 
instrument command and status words and editing 
the converted elevation angle. The conversion of a 
raw elevation value is dependent on the specific scan 
mode of the instrument. If the scan mode cannot be 
determined from the command and status words, or 
is changing during a record, no values of elevation 
angle are considered valid for the record; thus, no 
detector pointing vectors are computed. If the scan- 
ner is in the stowed mode, Earth viewing data are 
not available and no elevation angles are computed. 
Finally, before a converted elevation angle is used in 
computing a pointing vector, its value is tested and 
must be within specified bounds. 

Adjustment to detector elevation angle before 
computing pointing vector . Before the validated 

scan elevation angle can be used in the detector 
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pointing vector calculations, it is adjusted to account 
for scanner response time. In ERBE instrument 
studies at LaRC, it was concluded that the ERBE 
scanner radiometric detectors have a response time of 
46 milliseconds. This response time is greater than 
the 335-millisecond interval between consecutive ra- 
diometric measurements during the Earth viewing 
portion of a scan cycle. To account for this time 
lag, the scanner measurement locations are computed 
to match a set of radiometric detector outputs with 
the measurement locations that the detectors were 
viewing 46 milliseconds earlier. 

The four dynamic factors which affect the scanner 
Earth measurement locations are as follows: (1) ro- 
tational motion of the Earth, (2) spacecraft orbit 
motion, (3) spacecraft attitude motion, and (4) scan 
head elevation motion. At the equator, the Earth’s 
motion in 46 milliseconds is less than 25 m. The 
along-track spacecraft distance traveled in 46 milli- 
seconds is only 317 m for ERBS and 300 m for the 
NOAA-9 spacecraft. Changes in the location at nadir 
in 46 milliseconds as the result of attitude motion are 
normally less than 25 m and do not exceed 50 m. The 
ERBE scan head elevation beam rotates at a rate of 
66.7° per second during the Earth viewing portion 
of the scan cycle and thus rotates 3.066° in 46 milli- 
seconds. Earth distances measured at nadir due to 
this angular change are 32 km and 46 km for the 
ERBS and NOAA-9 spacecraft, respectively. Be- 
cause the combined changes in measurement location 
due to Earth rotation (0.02 km), spacecraft orbit mo- 
tion (0.3 km), and attitude motion (0.05 km) are very 
small in 46 milliseconds compared with the change 
due to elevation beam rotation (46 km), only the scan 
motion is considered in correcting the measurement 
locations. The actual adjustment in measurement 
location is made by computing the detector pointing 
vector in the spacecraft coordinate system using an 
effective detector elevation angle 3.066° less than the 
measured elevation angle. The most noticeable ef- 
fect of this adjustment is that the measurements are 
shifted toward the side where the scan begins. 

Transformation From Spacecraft Coordinates 
to Local-Horizon Coordinates 

This section describes how the detector pointing 
vectors determined in a spacecraft coordinate system 
are transformed into a common local-horizon system 
centered on the spacecraft. This local-horizon sys- 
tem is illustrated in figure 3, which also shows the 
relationship between the F-, and Z- axes of the 
system and the physical quantities which define the 
system. The X-axis is down and normal to the local- 
horizon plane, the Y -axis lies along the negative com- 
ponent of the spacecraft velocity vector in the local- 


horizon tangent plane, and the Z-axis is parallel to 
the orbit angular-momentum vector. A description 
of the orbit ephemeris data used to compute the 
position of the origin and the orientation of this 
coordinate system is given in the next section. 

The first step in transforming the detector point- 
ing vectors into local-horizon coordinates is to com- 
pute the pointing vectors in the local-horizon system 
for the specific spacecraft. The components of the de- 
tector pointing vectors in a spacecraft coordinate sys- 
tem and its components in the spacecraft-dependent 
local-horizon system are related through the space- 
craft orientation or attitude angles. (See fig. 4.) This 
relationship is different for the two spacecraft, both 
in the designation of the axes and in the manner and 
order of the rotations of the axes. 

The transformation of the detector pointing vec- 
tors to the ERBS local-horizon axis system, given its 
components in the ERBS spacecraft system, is given 
as follows: 


where 


'Xh,e' 


■x E - 

Yh,e 

— Ge^e^e 

Y E 

■ Z H,E - 


.z E . 



' COS0 

0 sin 0 ‘ 

E e = 

0 

1 0 


„ — sin 0 

0 cos#. 


'1 0 

0 - 

F e = 

0 cos <j> 

— sm<j> 


.0 sin</> 

cos <f> . 


'cos^> - 

- sin ip 0 ' 

G E = 

sin^ 

cosip 0 


. 0 

0 1. 


( 5 ) 


The transformation of the detector pointing vec- 
tors to the NOAA local-horizon system, given its 
components in the NOAA spacecraft system, is given 
as follows: 


’Xh,n‘ 


X N - 

Yh,n 
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II 
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( 6 ) 
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where 


vectors for the NOAA spacecraft are given directly 
in the coordinate system of figure 3 as follows: 


' cos# sin# 0 

Ejv = -sin# cos# 0 

0 0 1 

cos <t> 0 — sin 0 

0 1 0 

sin <f> 0 cos <(> 

1 0 0 
Gjv = 0 cos 0 sin 0 

.0 — sin 0 cos 0 

Sixteen sets of attitude angles (roll, pitch, and 
yaw) are contained in a telemetry record (one set 
per second). Since the attitude- angle sampling rate 
is much slower than the radiometric sampling rate, 
the pointing vectors for each second of radiometric 
data are computed using the set of attitude angles 
whose time most closely matches the average time of 
the 1-second group of radiometric data. Operations 
status data indicating the operational modes of the 
spacecraft are checked for every 16-second record, 
and the detector pointing vectors are computed only 
when the operational mode is explicitly determined. 
The attitude angles must be edited before they are 
used in computing the detector pointing vectors. 

The measurement location calculations are made 
for both spacecraft using data directly from the 
spacecraft onboard attitude determination systems. 
The Attitude Determination and Control System 
(ADACS) on the NOAA spacecraft is a proven sys- 
tem. It has been flown and verified on several differ- 
ent spacecraft. The attitude angles are usually quite 
small, and normally their values do not exceed 0.1° 
during a 24-hour period. The ADACS for the ERBS 
spacecraft is a new system, however, and the angles 
derived directly from the onboard system are larger 
than those for the NOAA spacecraft. Comparisons 
of attitude data from the ERBS onboard attitude 
system with data generated by more refined ground- 
based calculations indicate that the onboard system 
is sufficiently accurate to meet ERBE measurement 
location requirements. 

The local-horizon coordinate system in which 
the detector pointing vectors are calculated for the 
NOAA spacecraft by using the transformation ma- 
trix of equation (6) is identical to the common local- 
horizon system, in which all detector viewing vectors 
are finally determined. Thus, the detector pointing 



Xh i r xh,n ' 

y H = Y h ,n (7) 

Zh - 1 L Z H,N - 

The transformation equation (5) determines the 
detector pointing vectors in one of two ERBS local- 
horizon systems. The orientations of these two coor- 
dinate systems, which differ when the ERBS space- 
craft flies X-axis forward or X-axis rearward, are 
shown in relation to the common (NOAA) local- 
horizon coordinate system in figure 5. The ERBS 
spacecraft changes flight configurations by yawing 
(rotating about its Z-axis) 180° when the Sun is in 
the spacecraft orbit plane, which is the case approx- 
imately every 36 days. The transformations which 
determine the detector pointing vectors on the ERBS 
spacecraft in the common local-horizon coordinates 
are as follows: 

For X-axis forward 

'X H I r z HjE I 0 0 in rX HE - 

Yh = -Xjj,e = -1 o 0 Yh,E (8) 

-Zfj\ ~Yh,E J 0 -1 oJ Zfj,E- 

For X-axis rearward 

[Zh,e1 I" 0 0 il \ X H,E' 

Yh = X H,E = 10 0 Y HtE (9) 
-Zh- L l0 1 0J L Z H , E . 

Orbit Ephemeris Data 

The orbit ephemeris data for the ERBS and 
NOAA spacecraft are generated by GSFC. Tracking 
data for the ERBS spacecraft are obtained by the 
NASA Tracking and Data Relay Satellite and ground 
radar stations, and data for the NOAA spacecraft 
are provided by the U.S. Air Force North American 
Aerospace Defense Command. Eight days (one week 
with one overlap day) of ephemeris data are written 
on a single magnetic tape, which is shipped to LaRC 
weekly. Spacecraft position and velocity are con- 
tained on the ephemeris tape at 1-minute intervals, 
with the first data set for each day always at mid- 
night universal time. Each 8-day ephemeris data set 
is processed through software at LaRC which eval- 
uates the computed data for consistency and accu- 
racy before the data are merged with the correspond- 
ing telemetry data and used to compute scanner 
measurement locations. 
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Measurement Location Computation in 
Earth Fixed Coordinates 

In this second major step in determining the scan- 
ner Earth measurement locations, the unit pointing 
vector which describes the direction of the detector 
optical axis in the local-horizon coordinate system is 
transformed into an Earth fixed coordinate system. 
It is first necessary to transform all the pointing vec- 
tors and ephemeris data to the same coordinate sys- 
tem at the same time. This step is accomplished 
by merging instrument telemetry data and orbit 
position data as a function of time, interpolating the 
ephemeris data to the times of the scanner measure- 
ments, transforming pointing vectors and ephemeris 
states to the Earth fixed system, and finding the 
intersections of the pointing vectors with the Earth 
ellipsoid. 

Interpolation of Ephemeris Data Times to 
Telemetry Data Times 

As discussed previously, ephemeris state vectors 
are calculated at 60-second intervals, and telemetry 
data are recorded every 16 seconds. It is there- 
fore necessary to interpolate the ephemeris data to 
the time of the telemetry records. This step is ac- 
complished by converting two consecutive ephemeris 
state vectors to Keplerian orbit elements (ref. 2) and 
using these elements in linear interpolation to obtain 
orbit elements at the telemetry record times. These 
elements are then converted back to Cartesian state 
vectors to define the spacecraft position and velocity 
at the time of the telemetry record. Orbit element 
interpolation is justified because the Keplerian ele- 
ments, except for true anomaly, change very slowly 
with time. True anomaly varies almost linearly with 
time in the near-circular orbits used by the ERBS 
and NOAA satellites. 

Since scanner measurements are made at different 
times within a telemetry record, the spacecraft posi- 
tion must be further interpolated to the precise time 
of the measurement. The measurement time is ob- 
tained by adding the telemetry time offsets discussed 
previously to the universal time of the data record. 
Once the measurement times are known, the same 
linear interpolation scheme is used between teleme- 
try record times to find the spacecraft state vector at 
each measurement time. 

Coordinate Transformations 

Inertial coordinates to Earth fixed coordinates. The 
ephemeris data are given in the inertial coordinate 
system (fig. 6) with the X/-axis pointing toward the 
mean vernal equinox, the Zj-axis coincident with 
the Earth’s rotation axis, and the Yj-axis forming 
a right-hand system. 


The Earth fixed coordinate system has the 
Xjr-axis at the Greenwich meridian, the Z/r-axis 
along the Earth’s rotation axis, and the Yp - axis 
forming a right-hand system. The matrix for the 
transformation from inertial to Earth fixed coordi- 
nates is 


■x F - 


cos u sini/ O' 


xr 

Y F 

= 

— sin v cos v 0 


Yi 

. Zjr . 


1 

t-H 

O 

o 
1 


-Z T . 


where i/, the inertial longitude of Greenwich at the 
measurement time, is given in reference 2 by 

v = 99.6909833° + 36000. 7689°T U + 0.00038708°!^ 
+ 360.9856463° AT (11) 


where 


T u 

= (JD - 2 415 020) /36 525 = 
Julian centuries since 1900 

JD 

= Julian date 

A T 

= Fractional days between 
measurement time and 
midnight universal time 


Local-horizon coordinates to Earth fixed coordi- 
nates. Pointing vectors for both spacecraft are given 
in the local-horizon system shown in figure 3. The 
transformation from the local-horizon system to the 
Earth fixed system (fig. 7) is done in two steps. The 
first step transforms the pointing vector from the 
local-horizon system to the local-geodetic system, 
which has Xq pointing north, Yq pointing east, and 
Zq pointing “down.” The transformation matrix is 


x G ~ 


"0 —cos r] — sin?7 _ 


Xjj' 

Yg 

= 

0 sin 7] — cos r\ 


y h 

■ Zq. 


i—i 

o 

o 


,z H . 


where 77, the heading angle, is the angle between 
north and the projection of the velocity vector in the 
tangent plane. 

The second step of the transformation is from 
local-geodetic to Earth fixed coordinates as follows: 


1 

- 

Lz f J 

L 


— cos A sin 

— sin A sin 

cos 




— sin A 
cos A 
0 


— cos A$q 
— sin A cos $£ 
— sin $g 



'Xq- 


Yg 


-Z G . 

(13) 


8 


where A is longitude and is geodetic latitude. 
The relationship between geodetic and geocentric 
latitudes is found from the following equation (ref. 3): 

$G = tan -1 (tan $ c R e/R%) ( 14 ) 


where is geocentric latitude, is geodetic lat- 
itude, Re is the equatorial radius of the Earth, and 
Rp is the polar radius. 

The two transformations may be combined as 
follows to give the total transformation from local- 
horizon to Earth fixed coordinates: 


Xp " 


■ — cos A cos <f>Q cos A sin (j>Q cos 77 — sin A sin 77 cos A sin 4>q sin 77 + sin A cos 77 " 


-x H - 

Y F 

— 

— sin A cos sin A sin 4>q cos 77 + cos A sin 77 sin A sin <\>q sin 77 — cos A cos 77 


y h 

Zjr . 


— sin <f>Q — cos <f>Q cos rj — cos <f>Q sin 77 


■ z h ■ 


Location Geometry at Top of Atmosphere 

Once the spacecraft location and instrument 
pointing vectors are known in the Earth fixed co- 
ordinate system, the final step in locating the mea- 
surements is to find the intersection of the pointing 
vectors with the top of the Earth’s atmosphere. The 
Earth is modelled by an ellipsoid with polar radius 
Rp and equatorial radius Re . The height of the top 
of the atmosphere h is an arbitrary reference altitude 
(30 km) above the surface of the Earth and follows 
the equation 

X 2 /(R e + h ) 2 + Y 2 /{R iS + h ) 2 + Z 2 /(Rp + h) 2 = l 

(16) 

Figure 8 illustrates the geometry of the pointing 
vector intersection. Given a vector pointing to the 

spacecraft P and a unit pointing vector V, the 
magnitude of the vector from the spacecraft to the 

intersection with the ellipsoid V is found from a 
root of the quadratic equation with the following 
coefficients: 

A = (R P + h) 2 (V 2 + V 2 ) + (R e + h) 2 V 2 1 

B = 2 [(R p + h) 2 (P x \ x + PyVy) + (R e + h) 2 P z \z } I ^ 
C = (R p + h) 2 (P 2 +P 2 ) + (R e +h) 2 P 2 J 

- (Rp + h) 2 (Re + h) 2 

If B 2 — 4 AC < 0, there is no intersection. One of the 
real roots is given by 

I V| = ( -B - \ZW^1AC)/2A 

If this root is positive, it defines the magnitude of 
the vector from the spacecraft to the ellipsoid. If 


this root is negative, the pointing vector points away 
from the ellipsoid. The other root of the quadratic 
equation is the intersection of the pointing vector 
with the distant side of the ellipsoid. 

The vector R = P + V is the vector from the 
center of the Earth to the intersection point shown 
in figure 8. Finally, the latitude and longitude of the 
Earth located measurement in Earth fixed geocentric 
coordinates are 

= sin -1 (#«/! R|) 

A g = tan ~ 1 (R y /R x ) 

Method of Assessing Measurement 
Location Accuracy 

The technique presented in the preceding section 
is used to calculate the location of all radiometric 
measurements. Some method is required to assess 
the accuracy of the technique. One way to com- 
pare the accuracy is to identify Earth features in 
the radiometric measurements. The shortwave ra- 
diometric measurements are not used because they 
include large variations due to clouds and are not 
useful at night. When there are large thermal con- 
trasts, the longwave measurements show features 
such as a desert adjacent to an ocean. Using multiple 
scans across the thermal scene, the coastline between 
desert and ocean may be identified. This coastline 
may then be compared with an actual map, and er- 
rors in latitude and longitude or errors along and 
perpendicular to the subsatellite trace may be esti- 
mated. Coastline identification results for many days 
and several sites are accumulated on scatter plots, 
which allow biases in measurement locations to be 
detected. 
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Site-Selection Considerations and 
Restrictions 

Four measurement sites have been chosen to 
satisfy the following conditions: 

(a) Probable high thermal contrast between land 
and water. For example, at night during the summer, 
land is cooler than water, and during the day, land 
is warmer than water. 

(b) Infrequent cloud cover. When clouds are 
present in a coastline scene, the longwave scanner 
measurements are much colder than either land or 
water. These measurements can be removed from the 
data set by using a minimum temperature cutoff. 

(c) No unusual terrain features, such as lakes 
and mountains, next to the coastline. Lakes near 
the coastline have similar effects to those seen with 
ocean next to coastline, which makes discrimination 
difficult. Mountains may appear cold on one side and 
warm on the other side. This thermal contrast can 
be confused with that of land and water. 

(d) Interesting coastline. An absolutely straight 
coastline would not be useful in detecting errors 
along the coastline. A coastline with regular curves, 
peninsulas, and bays is most useful. 

In light of the preceding considerations, the 
following measurement sites have been chosen: 

(a) Baja, California 

(b) The northwest coast of Australia 

(c) The coast of Libya 


Lower measurements are assumed to be the reverse. 
Thus, if there is a homogeneous land region bounded 
by water, a scan across the coastline boundary yields 
measurements similar to those shown in figure 9. 

A cubic equation y * = ax ? + bx ? + cx{ + d (i = 
1,...,4) is fitted to each set of 4 consecutive points 
for latitude and longitude. The inflection point for 
each equation x = — b/Sa is calculated for each 
set of 4 points. These inflection points are consid- 
ered a coastal crossing if they fall between the two 
center points, X 2 and £3, and if the change in the 
measurement 

Ay = \yi ~ 2/4 1 

is larger than a threshold value taken here to be 
1 w/m 2 /sr. A threshold value is used to elim- 
inate spurious inflection points caused by normal 
variations in the measurements. 

The coefficients of the cubic equation are found 
from 


'a' 


t 3 

X 1 

T 2 

X 1 

XI 

r 

-1 

'yr 

b 


A 

x\ 

X2 

1 


V2 

c 


x 3 

x 3 

X3 

1 


V3 

.d. 


1 

H 

4^ Go 

x\ 

X^ 

1 . 


-y4- 


Once a point has been selected as a coastal 
crossing, its location is determined by interpolat- 
ing between the latitude and longitude of adjacent 
measurement locations. 


(d) The southeast coast of the Arabian peninsula 

Some restrictions are applied to the data. Data 
values less than a threshold are presumed to be 
clouds and are not used. To avoid atmospheric ef- 
fects such as limb darkening and refraction, data 
with scanner pointing angles greater than 30° mea- 
sured from nadir are eliminated. To reduce extra- 
neous thermal contrast due to inland terrain, data 
estimated to be more than 25 km from the coastline 
are not considered. 

Finding a Scanner Coastline Crossing 

The intersections of the pointing vectors with the 
TOA are extrapolated to the surface of the Earth 
(fig. 8), and these points on the surface are examined 
for coastline crossings. Figure 9 shows the plot of 
a typical longwave scan as it crosses a coastline. 
Each point represents one longwave scanner radiance 
measurement. Higher thermal measurements are 
presumed to be land in daytime and water at night. 


Fitting a Collection of Coastal Crossings to 
a Map 

Once a collection of coastal crossings has been 
calculated, it is fitted to a digitized map of the 
coastline of interest. It should be noted that the map 
coordinates are geodetic coordinates which differ in 
latitude from the calculated geocentric coordinates. 
This relationship between the two latitudes is given 
in equation (14). 

Figures 10(a) through 10(d) show good examples 
of sets of crossings plotted on maps. The circles 
represent the calculated coastline crossings, and the 
squares are the subsatellite points. Most sets of 
crossings do not produce useful information because 
of the presence of clouds or poor thermal contrast at 
the time of satellite passage. 

Some figure of merit is needed to determine how 
well the estimated coastline matches the map coast- 
line. The figure of merit chosen is the least-squares 
distance between the calculated coastal crossings and 
the map coordinates of the coastline. A latitude and 
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longitude correction is applied to each member of 
the set of crossings in a manner which minimizes the 
sum of the squares of the distances from the crossings 
to the map. This correction is then defined as the 
location error for the satellite pass. 


Along-Track and Cross-Track Errors 

Since much of the error in locating measurements 
may be attributed to spacecraft attitude errors and 
errors in the adjustment to detector elevation angle, 
it is useful to consider the measurement location 
errors in a coordinate system which is aligned with 
the scan direction. The mapping from latitude and 
longitude coordinates to along-track and cross-track 
coordinates is shown in figure 11 and is given by 


£c" 


± sin 7} cos r/ 

ex 

^a_ 


— cos rj ± sin 7] 

, € <s>. 


( 20 ) 


where r/ is spacecraft heading angle, e c is cross-track 
error, e a is along-track error, is longitude error, 
and €$ is latitude error. The upper sign (+) in equa- 
tion (20) is used when the scan direction is left to 
right; otherwise, the lower sign (— ) is used. The lati- 
tude effect on longitude error is ignored, because the 
sites considered are near the equator. For example, 
figure 10(a) shows that the least-squares latitude and 
longitude errors map into errors of —0.0478° cross 
track and 0.0271° along track. 


Statistical Treatment of Results 

The along-track and cross-track errors from dif- 
ferent sites and different days may be used to de- 
termine biases in the location of the measurements. 
Figure 12 shows a scatter plot of data from the ERBS 
satellite for November 1984. Each point on the plot 
is the result of analyzing the coastline crossings for a 
single site. The solid ellipse is the 95-percent proba- 
bility ellipse computed assuming that the data have 
a Gaussian distribution with mean and covariance 
equal to that of the data. The center of the el- 
lipse identifies a bias along each axis. Figure 12 
shows a cross-track bias of —0.026° (about 2.9 km) 
and an along-track bias of 0.0048° (about 0.53 km). 
The variance implied by the probability ellipse is an 
underestimate, because scan angles are limited to 
those within 30° of nadir. For a given pointing error, 
the error in Earth location increases with scan angle 
as measured from nadir. 


Results and Discussion 

The method for calculating ERBE scanner mea- 
surement locations presented here has been imple- 
mented in the operational data processing system. 


Every scanner measurement is located by using the 
full set of transformations and interpolations. These 
results are based upon operational flight data. 

The coastline detection method has been applied 
to many days of data for both the ERBS and NOAA 9 
satellites. Figures 10(a) through 10(d) illustrate 
orbital passes over Baja, California, the northwestern 
coast of Australia, the Gulf of Sidra off Libya, and 
the southeast coast of the Arabian peninsula. The 
square symbols show the subsatellite track during 
the pass; the detected coastline crossings for each 
scan are drawn as circles. The estimated coastlines 
appear to follow the map quite well. In figure 10(a), 
one point in the northeast portion of the Gulf of 
California appears to be 20 to 30 km from the coast. 
In fact, there is an island, Isla del Tiburon, which 
is not digitized on the map but is detected by the 
scanner. 

Individual scanner measurements in figure 10(a) 
show a difference of up to 15 km between the es- 
timated coastline crossing and the actual mapped 
coastline. When the coastline fitting technique is 
used, the average error for all crossings during the 
pass is 6.1 km. This error may be apportioned be- 
tween 0.027° along-track error and —0.048° cross- 
track error. These errors are approximately 3 and 
5 km, respectively. 

The location error for a single satellite pass over a 
coastline may not be representative of the operation 
of the satellite over all orbital conditions. Therefore, 
additional passes must be investigated. A scatter 
diagram of along-track and cross-track errors for all 
coastline crossings identified for the ERBS satellite 
during November 1984 is presented in figure 12. Each 
point in this diagram represents the average error for 
a single map of one of the coastal areas. For example, 
the results for figure 10(a) are collapsed to a single 
point in this plot. A small but persistent cross-track 
bias is illustrated by the 95-percent confidence ellipse 
which is centered about —0.03° in the cross-track 
direction. This direction is opposite to the direction 
of scan motion during this time period. The bias 
represents a 3-km location error. To correct for this 
bias, the instrument response time adjustment to the 
effective scan elevation angle has been reduced from 
3.066° to 2.85°. 

Figure 13 shows results for the northwest coast of 
Australia. Figure 13(a) uses a response time adjust- 
ment of 3.066°, and figure 13(b) uses an adjustment 
of 2.85°. A comparison of the figures shows that 
the estimated coastline is closer to the map coast- 
line for the reduced instrument response time adjust- 
ment. The corrected time adjustment is used in the 
operational processing system. 
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As discussed previously, the ERBS Attitude 
Determination and Control System is new, and 
the attitude angles derived from the onboard sys- 
tem are larger than those for the NO A A space- 
craft. Ground-based processing of ERBS AD ACS 
data has been performed at GSFC to produce more 
accurate (definitive) attitude data for selected pe- 
riods. The coastline detection results for ERBS 
shown in figure 14 illustrate scanner measurement 
locations based on a set of these definitive attitude 
data. The results shown are for the same geographi- 
cal area and time period as those of figure 10(d), but 
the scanner measurement locations in figure 10(d) 
are based on attitude data derived from the onboard 
attitude determination system. This comparison is 
typical, and the very close agreement between the 
two sets of coastline detection results indicates that 
the attitude data from the ERBS onboard deter- 
mination system is sufficiently accurate for ERBE 
measurement location calculations. 

Figure 15 presents April 1985 results for both 
the ERBS and NOAA 9 satellites. Again, the dia- 
gram plots along-track location errors versus cross- 
track location errors by using the average error dur- 
ing a satellite pass. These errors include the bias 
correction to the response-time adjustment for both 
satellites. The error ellipsoids show that 95 percent 
of all coastline detections fall within approximately 
10 km of the true coastline location. The mean and 
standard deviation along each axis of these scatter 
diagrams are shown in the following table: 



Cross track 

Along track 



Standard 


Standard 


Mean, 

deviation, 

Mean, 

deviation, 

Satellite 

km 

km 

km 

km 

ERBS 

0.09 

4.61 

0.49 

3.67 

NOAA 9 

-.39 

4.88 

-1.17 

3.43 


The average errors are much smaller than the 
scanner field of view, which is about 40 km. Also, the 
scatter of points as measured by the standard devia- 
tion along each axis is consistent between satellites. 
This consistency shows that error sources contribut- 
ing to the measurement location errors are similar 
for ERBS and NOAA 9. 

Conclusions 

Based upon the results from the coastline detec- 
tion technique developed here, the accuracy of the 


Earth Radiation Budget Experiment (ERBE) scan- 
ner measurement location appears to be well within 
the limits of the detector field of view. Results of this 
study indicate the following conclusions: 

1. Satellite ephemeris data are accurate and are 
merged properly with telemetry data for both the 
Earth Radiation Budget Satellite (ERBS) and the 
NOAA 9 satellite. 

2. The ERBS onboard attitude determination 
system produces attitude angle data which are suffi- 
ciently accurate for ERBS scanner location calcula- 
tions. 

3. At pointing angles within 30° of nadir, single 
estimates of a coastline location from the scanner 
radiometric data are less than 15 km in error. 

4. An average of all coastline detections from a 
single satellite pass is less than 10 km in error. 

5. For both the ERBS and NOAA 9 satellites, all 
passes of April 1985 have an average location error of 
less than 1.2 km, and the standard deviation of the 
error is less than 5 km in both the along-track and 
cross-track directions. 

Users of the ERBE data may rely upon the 
geographic location as identified on the archival prod- 
ucts. However, it should be noted that the geocen- 
tric coordinates identified in the ERBE data-tape 
products must be transformed to geodetic map co- 
ordinates before making comparisons. The ERBE 
project intends to monitor the measurement lo- 
cation accuracy by using the coastline detection 
technique for the life of the instruments. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
June 18, 1987 
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(a) ERBS. 



(b) NOAA. 

Figure 4. Relationship between a spacecraft coordinate system and corresponding local-horizon coordinate 





Figure 5. Relationship between ERBS local-horizon coordinate systems and final local-horizon system. 
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Figure 8. Intersection of pointing vector with Earth ellipsoid. 
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Figure 9. Cubic fit for coastline detection. 





Figure 10. Continued. 






Figure 10. Concluded. 
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Figure 13. Effect of instrument response time adjustment to effective scan elevation angle for Australian site. 
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Figure 13. Concluded. 







(a) ERBS. 

Figure 15. Ensemble location errors for April 1985. 
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(b) NOAA 9. 
Figure 15. Concluded. 
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